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The baryon density which may be produced during the electroweak phase transition in 
supersymmetric models is computed, taking into account the previously neglected effects 
of transport, strong and weak anomalous fermion number violation, thermal scattering, and 
a new method for computing CP violating processes during the transition. We can account 
for the observed baryon asymmetry, provided new CP-violating phases are greater than 
~ 10 - ( 2-4 ), and some superpartners are light enough to be relevant during the transition, 
which takes place at a temperature of (50-100) GeV. In one case, light superpartners are 
the top squarks and the charginos and/or the neutralinos; in another case the top squarks 
and both Higgs doublets are light. Our calculation is easily extended to the case of a 
general two Higgs model, where we find sufficient baryogenesis provided that a certain 
combination of parameters in the Higgs potential leads to a CP violating space dependent 
phase in the top quark mass of order 10~ 3 . 
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1. Introduction 



Supersymmetry is an attractive candidate for the physics of electroweak symmetry 
breaking, while electroweak baryogenesis (EWB) 0] is an explanation of the origin of 
the cosmological asymmetry between matter and antimatter in terms of experimentally 
accessible physics. It is therefore of interest to understand whether EWB is feasible in 
supersymmetric models. The only previous estimates []2],[| of the baryon asymmetry pro- 
duced in supersymmetric models neglected many effects which are now understood to be 
important, such as transport [f§-|| and thermal scattering [0-|9[|. 

Let us review the physics relevant for EWB. Anomalous baryon violation in the weak 
interactions takes place via unobservably slow tunnelling processes at zero temperature 
ICQ , but at temperatures above the critical temperature for the weak phase transition, 



theoretical estimates give a rate T = na^T, where a w is the weak fine structure constant 
O], and k is a pure number of order one@. Thus electroweak baryon number violation 
is fast enough in the early universe to change the cosmological baryon number. In ther- 
mal equilibrium, unless some nonanomalous approximately conserved quantum number is 
nonzero ||14|1 , anomalous processes will wash out any net baryon number, however a first 



order electroweak phase transition can provide the departure from thermal equilibrium 
necessary to generate a nonzero baryon number. Electroweak baryogenesis is only feasi- 
ble if two conditions are met, which probably require new weak scale physics beyond the 
Minimal Standard Model (MSM) ||J|, 15] - 19[ . (For relatively recent reviews, see [|0)-) 



1. The transition must be strongly enough first order so that after the transition the 
anomalous baryon number violation is too slow to wash out the baryons created during 
the transition ||21[| . This rate is proportional to exp(— M s /T), where M s , the energy 
of the sphaleron field configuration, is proportional to the W boson mass, M s = 



(90 — 160)M W ||22|1 . The condition that the W mass jumps to a large enough value 
during the transition to avoid post-transition baryon number washout requires a light 
Higgs in the MSM ]T8|j2^ , pi| (in lattice simulations, the transition appears too weakly 



1 An early estimate of k gave k > 0.1 [12] and a recent computation claims a value of k 



1.09 ±0.04 |3|. 



first order unless m# <C M w [^5] ). However with a top mass of (170-200) GeV, if 
the MSM is valid up to 10 6 GeV we will only be in the MSM ground state today for 



a Higgs mass heavier than ~ M w [26 



2. The amount of CP violation must be just right to explain the observed baryon to 
entropy ratio, ris/s ~ 10 -10 . The CP violation in the minimal standard model is 
only physical in processes which involve all the Cabbibo-Kobayashi-Maskawa (CKM) 
angles and in which all the like charge quark mass differences play a role, which makes 
it seem a priori difficult for Kobayashi-Maskawa CP violation to generate sufficient 
baryon number during the weak transition. An interesting attempt to find a large en- 
hancement of the CKM contribution to electroweak baryogenesis was made by Farrar 
and Shaposhnikov [^7| but was later shown not to work due to quantum decoherence 
effects 0,| . 

In constrast to the MSM case, in most extensions of the standard model there can be 
additional sources of CP violation which appear in particle mass matrices. During a first 
order electroweak phase transition, bubbles of the broken phase nucleate and expand. In- 
side the bubble wall, particle mass matrices acquire nontrivial space-time dependence and 
cannot be made real and diagonal at all points without introducing new CP violating terms 
into the particle dispersion relations. In a recent paper || we introduced a general method 
for computing the effects of the CP violating mass terms on particle distributions, which 
takes into account both the effects of scattering from thermal particles and the terms which 
lead to CP violation in particle propagation. It is now established that transport of CP vi- 
olating quantum numbers into the symmetric phase, where anomalous electroweak baryon 
number violation is relatively rapid, plays a dominant role in electroweak baryogenesis for 
all bubble wall widths i-iH|JT|]. 

The most well motivated viable theories for weak scale baryogenesis are two (or more) 



Higgs models[]28[ and models with weak scale supersymmetry. In the two Higgs model the 
relevant CP violation is produced by a phase in the Higgs potential, which leads to CP 
violating mass matrices for fermions and Higgs bosons, and produces especially large CP 
violating effects on the Higgs and axial top number distributions. Experimental constraints 
on atomic and neutron dipole moments allow the relevant phase to be as large as 0(1) 



|29fl . Also, the two Higgs model can easily simultaneously satisfy the constraints on Higgs 



particle masses and the requirement of a sufficiently first order transition []28| , j30l , |3~l| . There 
are many possible supersymmetric extensions of the MSM, with additional CP-violating 
phases. The minimal additional particle content (the Minimal Supersymmetric Standard 
Model, or MSSM) includes superpartners for all particles and a second Higgs doublet. The 
supersymmetric terms in the Lagrangian do not introduce any additional CP violation, 
however super symmetry must be broken by adding soft supersymmetry breaking operators, 
which in general are CP violating. If the new CP violating phases are of order one, 
the neutron and atomic electric dipole moments are larger than the experimental bounds 
pm unless the superpartners are unnaturally heavy [|33]], hence the usual assumption is 
that the soft supersymmetry breaking terms arise from CP conserving physics and have 
negligible phases. However it has recently been argued that in most grand unified 
supersymmetric theories, renormalization of the soft operators between the Planck mass 
and the scale of grand unification will induce phases of order 10 -2 — 10 -3 in the soft susy 
breaking operators, providing a new source of CP-violation into the low energy effective 
theory, which is just beyond our current experimental reach. In this paper we will assume 
that the supersymmetry breaking terms have CP-violating phases and see whether these 
phases can account for sufficient baryogenesis without violating the electric dipole moment 



bounds [|35[ . In the MSSM, the requirement of a sufficiently first order phase transition 
places upper limits on the Higgs and stop (supersymmetric partners of the top) masses 
||36|| which are barely consistent with experimental constraints-there are speculations that 
these bounds could be relaxed slightly by higher order and nonperturbative effects [|37|P. 
The MSSM may easily be extended by adding a gauge singlet |3£| which substantially 
removes these constraints. Here we will consider models both with and without a singlet, 
but we will only consider those sources of additional CP violation which may be present 
in the MSSM, with a general set of soft supersymmetry breaking terms consistent with 



experimental bounds. Therefore we will not worry about the mass upper bounds of ref. |36 



2 Note that in supersymmetric models, the vacuum stability lower bounds p6| ] on the Higgs 
mass do not apply. 
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and we will assume the Higgs potential is CP-conservingfl We refer to the supersymmetric 
model either with or without additional gauge singlets as the Supersymmetric Standard 
Model, or SSM. 

In the next section we discuss the dominant baryogenesis mechanisms in the SSM. 
In §3 of this paper we write down the set of coupled differential equations which describe 
particle interactions and transport during the weak phase transition, and make reasonable 
approximations which allow us to find an analytic solution for the baryon asymmetry in 
the SSM. In §4 we do the same for the two Higgs model. We conclude with a summary of 
our results and their implications in §5. 



2. CP violation and Particle Sources in the SSM 



Following previous work [|]-|||T5|-|T9| we compute the baryon asymmetry using the 
following steps: 

I. Compute the CP-violating perturbations of the plasma locally induced by the passage 
of the wall ( "particle source terms" ) . In ref . |§ we described all the sources in terms 
of quantum mechanical CP-violating reflection and transmission from layers of the 
phase boundary, combined with re-thermalization of the phase-space distributions. 
Unlike earlier calculations, whose applicability was restricted to either a "thin wall" 
or a "thick wall" , referring to whether the wall thickness is larger or smaller than the 
relevant mean free paths, our approach provides a unified and consistent treatment for 
all values of the wall thickness. The proposed method links the charge generation to 
microphysical CP violating processes, and hence can be widely applied. It generalizes 
the method developed in Ref. |7J and so properly incorporates decoherence effects 
which have been shown to have a major negative impact on the generation of a CP 
violating observable in the MSM [7||]]- 
II. We approximate the solution to the Boltzmann equations for particle distribution 
functions by writing down and solving a set of coupled differential equations for the 



3 In some models with a singlet there can be CP-violation in the Higgs potential which can 
produce CP-violating effects very similar to those in two Higgs models, however in most models 
the Higgs potential automatically conserves CP. 



local particle densities including the source terms, transport, Debye screening [^9 



of induced gauge charges^ and particle number changing reactions ||. The solution 
to these equations generally includes a net baryon number, which is produced in the 
symmetric phase and is transported into the bubbles of broken phase, where it survives 
until the present provided that the phase transition is sufficiently first order. 
In this section we focus on the first step in the calculation. CP-violating particle 
source terms have been shown to result for a selected subset of species in the plasma which 
mix with one another via a mass matrix with complex phases which either: 



a) cannot be rotated away as the result of interactions with the plasma [27[ . 

b) cannot be rotated away at two adjacent points x and x + dx, by the same set of unitary 
transformations, that is, U~ U x +dx 7^ 1- 

When present, the second mechanism dominates over the first one, as the first mecha- 
nism generically involves additional particles whose coupling to the plasma, yields further 
suppressions. It is the second mechanism which controls baryon generation in the SSM as 
the neutralinos, charginos, and squarks, have mass matrices with CP violating entries and 
a non-trivial space-dependence due to the Higgs vacuum expectations values (c/.§2.2). So 
is the case in the two Higgs models with explicit CP violation in the Higgs potential which 
yields a space dependent phase to the top quark and Higgs massesi (c/.§4). In contrast, 
in the minimal standard model, the quark mass matrix has only an overall dependence on 
the Higgs vacuum expectation value and can be diagonalized by space-independent unitary 
rotations, hence, it can generate a CP violating observable through mechanism a) rather 
than through mechanism b) - z.e.through charge current interactions which correct the 
dispersion relation of the propagating quark in the plasma ||27|1 . This mechanism, however, 



has been shown to be quite ineffective at generating a significant baryon asymmetry in the 
MSM 0,1. 



4 In practice we simplify our equations by ignoring the effects of screening since the impact on 
baryogenesis turns out to be of order 1 [40,41]. 



5 Another potentially relevant species in this model is the r-lepton, and some have argued 
that its contribution dominates that of the top quark and Higgs |31,4^|. We do not confirm the 
importance of the r lepton unless tan (3 is very large. 



The method introduced in Ref. can account for both mechanisms. However, as we 
are concerned with extensions of the standard model for which the second mechanism b) is 
dominant, we will review the general principles of the method for this specific situation!. 



2.1. The Method 

Let us consider a set of particles with mass matrix M(z) and moving, in the rest frame 
of the wall, with energy-momentum E, k. We wish to find the CP-violating asymmetry 
in their distributions which results from their passage across the wall. We define z D to be 
their last scattering point, where they emerge from a thermal ensemble with a probability 
distribution represented by a density matrix p Zo . These particles propagate freely during 
a mean free time r, then rescatter and return to the local thermal ensemble in the plane 
z a + tv, v being the velocity perpendicular to the wall, k±/E. During the time r, these 
particles evolve according to a set of Klein-Gordon, Majorana or Dirac equations coupled 
through the mass matrix M(z). (Some effects of interaction with the plasma, which do not 
destroy quantum coherence, can easily be included in these equations.) In the course of this 
evolution CP violation affects the distribution of these particles. At z a , the contribution 
of these particles to any given charge cancel exactly the contribution of their antiparticles 
since the charge is CP odd and we take the density matrix to be CP even. However, after 
evolving a time r across the CP violating space-dependent background, this cancellation 
no longer takes place. At the subsequent scattering point z a + rv , these charges assume a 
non-zero value, as the evolution of the particles over the distance rv can be CP-violating. 
Specifically, the probability for a particle emitted at zq to be transmitted to zq + rv can 
be different from the transmission probability for its CP conjugate. It is only necessary 
to follow the contribution of a selected subset of charges carried by these particles in 
order to characterize completely the departure from thermal equilibrium resulting from 
the passage of the wall. In this subset, there are charges which are explicitly violated by 
the mass matrix (axial charge . . . ) and charges which are exactly conserved ( B-L, . . .). 
In the latter case, there is no net charge density created but, instead, a net spatial current 

6 Our method with mechanism a), applied to the standard model, would give results in agree- 
ment with the ones obtained in ref. JjJ, where similar techniques have been used. 
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emerges: opposite charges move in opposite directions. This spatial current arises at high 
order in masses and the wall velocity v w and, in most of the cases we considered, yield a 
subleading contribution to the baryon asymmetry. We are mostly interested in the former 
charges, that is, those charges which are violated by the presence of the wall. In addition 
to a spatial current, those charges develop a net average density; the latter is linear in the 
wall velocity and arises at low order in the mass expansion. We will focus our attention on 
those charges which develop a net density. We alert the reader that there may be situations 
in which conserved charges have to be equally considered. 

To be more quantitative, we introduce J±, the average current resulting from particles 
moving toward positive(negative) z between z a and z Q + A, where 

A = rv. (2.1) 

The current J_|_ receives contributions from either particles originating from the thermal 
ensemble at point z Q , moving with a positive velocity and being transmitted at z a + A, 
or from particles originating at z Q + A, moving with velocity — v and being reflected back 
towards z a + A (Fig. la). A similar definition exists for J_(Fig. lb). J± are CP violating 
currents which are associated with each layer of thickness A moving along with the wall; 
they can be computed according to 



J+ =! { Trp Zo [T^QT - T f QT] + Trp Zo+A [&QR - R QR] 



J_ 



R ] QR - R ] QR 



-t 



f ] QT - f QT 




(2.2) 



R(R) and T(T) are reflection and transmission matrices of particles (anti-particles) 
produced at z Q with a probability matrix p Zo , evolving toward positive z (increasing mass); 
R and T are the corresponding quantities for particles produced at z Q +A with probabilities 
contained in p Zo +A and evolving toward negative z; v is the magnitude of the group velocity 
perpendicular to the wall at point z Q while v is the same quantity but a distance A away. 



Finally, Q is the operator corresponding to the chosen charge and the trace is taken over 
all relevant degrees of freedom and averages over the location z a within a layer of thickness 
A. 

Formulae (|2.2|) provide a concise method of computing the CP violating charge cur- 
rents J±, which results from the propagation and the mixing of particles within the layer 
of thickness A at point z a . After a boost to the plasma frame, these currents constitute the 
fundamental CP violating building blocks that we need to construct the source terms of 
the system of rate equations introduced in §3, which ultimately will convert them through 
diffusion and relaxation mechanisms into a net baryon asymmetry. 

For our purpose, we construct the source terms as follows. Consider a small volume 
element in the plasma. As the wall crosses it, it deposits into it the current density 
( J_i_ + J-) pi asma every time interval r; the subscript ) p i aS ma refers to the quantity boosted 
to the plasma frame. At an arbitrary time t, the current density accumulated by this 
mechanism is@ 



f \dt'(j + ^t') + J.^t% lasma . (2.3) 

Jt-TR T 

Here, tr is a typical relaxation time. From this, we infer the net rate of change of charge 
Q per unit volume to be 



7 Q (x,f) =d^s 



— ~ (J+(%i t) + J-(x, t))p lasma — — (J+(x, t — Tr) + J-(x, t — TR))p lasma ^ 4) 



f 1 

/ ~dz(J+ + J—)plasma- 

Jt-T R T 



Formula ( |2.3|) along with formulae ( |2.2| ) constitute the starting point for our analysis 
of the SSM and of the two Higgs model. In practice, we simplify eq. ( [2.4| ) by making an 
expansion in spatial derivatives, allowing us to neglect the third term, and we take tr large, 
so that we can neglect the second term, whose effect is accounted for as an independent 
relaxation term in our rate equations (cf. §3). 



7 We leave aside diffusion which is accounted for independently in the rate equations. 



One particular advantage of the method above is that it does not require any assump- 
tion about the relative magnitude of the mean free paths and the thickness of the wall. 
Hence, in contrast with earlier methods, it unifies all electroweak baryogenesis scenarios. 



2.2. The SSM 

In the SSM, we are interested in the generation of charges which (a) are approximately 
conserved in the unbroken phase so that, they can diffuse a long way in front of the 
bubble wall, where anomalous baryon violation is fast and (b) are non-orthogonal to baryon 
number, so that their relaxation energetically favors a non-zero baryon charge. Candidates 
of choice are Higgs number and axial top number. The generation of these charges results 
from the mixing of the charginos, neutralinos and the mixing of top squarks respectively. 
The chargino mass matrix is, in the basis W + , W_, h_, h' + 



( 

-V 2 



m 2 





-v 2 





\ 

-vi 



(2.5) 



V -v x e^ B ii / 

with v x = y/2M w (z,T) cosp(z) and v 2 = V2M W (z, T) smp(z); M w (z,T) is the 
temperature-dependent W mass defined at each point z in the wall. The neutralinos 
mass matrix is, in the basis W3, B, ho, h' Q , 



Mn = 



( m 2 


u 2 cos 9 W 
V —u\ cos 9 U 





fhi 
—u 2 sin 9 W 
u\ sin 6 W 



u 2 cos 6 
—u 2 sin 9 




w 
w 



-u\ cos 9 W \ 
u\ sin 9 W 



-e l(pB ix 







(2.6) 



/ 



with u\ = Mz(z,T) cos j3{z) and u 2 = Mz(z,T) sin j3{z). Finally, in the basis t^, t^, 



ml 

a 2 e~ lct 



a 2 e ia 



m 



(2.7) 



here, we defined a 2 e ia = v^y {Ae i<t>A v 2 + Lie^ B v x ). 

With these conventions, the charge operator for the Higgs number takes the form 



Q~ h = Diag(0, 0, 1,-1), 
9 



(2.8) 



while for the axial stop number, the charge operator is defined to be 

Qs = Diag(I, -I). (2.9) 

We can now compute sources for those charges. We choose to perform an expansion in 
powers of mass. This expansion, introduced in @ and further developed in ||, is adequate 
to demonstrate the quantum mechanical physics required for generating a CP violating 
observable. In particular, it generates polynomials in M whose imaginary part of the 
trace yields an expansion in terms of CP violating invariant M. We will discuss the validity 
of this approximation later on in this section. 

- Charginos and Neutralinos - 

In order to compute the source 7^(x, t) for the Higgsino number, we begin by com- 
puting the corresponding current sources J± as given in Eqs. ( |2.2| ). Their determination 
requires the knowledge of transmission and reflection amplitudes which are obtained by 
solving a set of coupled Majorana equations with the chargino and neutralino mass matri- 
ces given in fl2.5| )-( |2~lf ). We obtain up to an overall phase, at leading order in M 2 ~ 

T =1 - [ dzi [ 1 dz 2 M 2 Mle i2u{zi - Z2) 
Jo Jo 

+ [ dz x [ 1 dz 2 [ dz 3 [ 3 dz±M 4 M* 3 M 2 M\ e *M*i-*a+*3-*4) + _ 

Jo Jo Jz 2 Jo 



(2.10) 



r A 

R = dzxMxe' 1 ^ 
Jo 

+ [ dz 2 [ dz 3 [ 3 dz 4 M 4 M* 3 M 2 e l2 -(-^3-. 4 ) + 

Jo Jz 2 Jo 



(2.11) 



and similar expressions for T, R. The quantity uj stands for the magnitude of the energy 
of motion transverse to the wall, and Mi is short for M(c,n){ z i)- The leading CP violating 
contributions arise at order 0(M/T) 4: . Contributing paths are depicted in Fig. 2. 

In order to use Eq. (|2.2| ) to compute J±, we need the density matrices p Zo and p Zo +A, 
describing the distribution of particles in phase and flavor spaces at their production point 



There are also CP violating self-energy corrections, which are the main source of CP violation 
in the absence of the one considered here. 
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z Q and z Q + A. We choose the density matrices describing equilibrium distributions of mass 
eigenstates in the unbroken phase 



P n Zo =~Dia-g(n^ 2 (E,v) , n^{E, v) , n^(E,v) , n^E, 



(2.12) 



We construct pI^+a f rom P^ n with the substitution v <-» — u . In these expressions, n m (E, v) 
is the Fermi-Dirac distribution, n/, of a species of mass m, boosted to the wall frame, 



71, 



(ExpfyviE - v w k ± )] +1) 1 • 



(2.13) 



Some motivations for our choice of density matrix are as follows. In a regime of large 
masses fh,2, p > T, it is obviously sensible to assume that p is diagonal in the mass 
eigenstate basis; in a regime of small masses, < T, particles are produced as interaction 
eigenstates which differ from mass eigenstates by a unitary rotation; ignoring this rotation 
amounts to ignoring small corrections of order (M/T) 2 . Furthermore, the choice of a 
thermal distribution is reasonable as the non-equilibrium component of the distributions 
is of order v w , that is, it amounts to ignoring terms of order < 1 pi,|I3]. Inserting Eqs. 
( ETT0D and ( FTTH) in Eq. Q yields 



+ 



(1, 0, 0,5) x | 



-A pzi pA pz 3 

+ 4 / dz\ / / dzs / dz^ sin 2uj (z4 — zs + Z2 — zi) 

J z 2 















xlm Tr 



p 2o M 2 Ml M A Q~ h 



"A pA pA pZ3 

— 4/ ctei / dz2 / ^3 / dz4 sin 2u (z4 — zs + Z2 — zi) 

J Z2 







xlm Tr 



o 



o 



Ml p Zo +A M 2 Ml M A Q- h 



[2.U) 
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and 



J_ = (1, 0, 0,-v) x | 

pA pA pz 2 pA 

— 4 / dz\ I dz2 / dz3 / dz4sin2uj(z4 — zs + Z2 — zi) 

JO Jz! JO Jz 3 

p Zo +A Ml M 2 M* 3 M 4 Q- h 

/•ii pA pz 2 pA 

+ 4 / dz\ I dz2 / dzs / dz4sin2uj(z4 — z% + Z2 — zi) 

Jo Jo Jo Jz 3 

M{ p Zo M 2 Ml M± Q~ h 



xlm Tr 

"A 



xlm Tr 



(2.15) 



In order to proceed with analytic expressions, we simplify further by performing a deriva- 
tive expansion: M(z) = M(z ) + (z — z Q ) d z M(z ) + 0(t/w) 2 and v = v + 0(t/w) 2 . 
This expansion is only justified in a region of the parameter space for which the mean 
free time r is smaller than the scale of variation of the masses, i.e., the wall thickness w. 
We will discuss the validity of our approximations in the last section. This simplification 
allows us to perform the trace in flavor space and in phase space independently and yields 
after summing the contributions from charginos and neutralinos, in first order in the wall 
velocity v w , 



(J+ + J_) 2 = 
(J+ + J_)° = lw v w x J* 



eigen. 



X 



d 3 k fjiUjA) r 

(2*)* {2Vi) 



E 



f T E t (2-16) 



(l + e Bi/T)2 T 



with /(£) defined as 



/(0 = (£cos£ - sinO' 



(2.17) 



and the summation is over the mass eigenstates with eigenvalues mi in the unbroken phase, 
with J" 1 , a corresponding CP-violating invariant. For the charginos, there are two mass 
eigenstates with masses \x and rh-2, whose CP-violating invariants are 



Jch =J2 2 = - Sin ^ V ™2 (V2d z v! ~ V!d z V 2 )/T 5 

= 2 n m 2 M$y(T, z) sin0 B d z (3 /T 5 . 
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(2.18) 



While, for the neutralinos, there are four eigenstates with masses fi, fj,, mi and m 2 and 
with, correspondingly, 

J™ 1 = - sin^s /U rhi sin 2 9 W {u 2 d z u\ - u 1 d z u 2 )/T 5 

= (irhi M 2 z (T,z)s\n 2 6 w sin0 B d z (3 /T 5 (2.19) 

J™ 2 = - sin S jx m 2 cos 2 9 w (u 2 d z ui - Uid z u 2 )/T 5 
= ixm 2 M Z (T, z) cos 2 d w sm(p B d z (3 /T 5 . 

We can now construct the local source density J^x, t) for the Higgs number by in- 
serting Eq. ( |2.16| ) in formula ( |2.4| ), we obtain 



rp4 

+0{v 2 w ,(t/w) 2 ) 

The factor Z~™ contains information on the phase space as well as on effects due to 
plasma interactions. Its analytic form is 



r oc y r zv-(9) a )/t, f[yt +t^T 

T£ = V^T / dy y 2 -r / d£ —± (2.21) 

h 7^ (l + e^) 2 J ^ y + m^f 2 K } 

it is a simple exercise to show that the factor Z~™ vanishes rapidly with r, as 

Jf oc r 5 (r->0). (2.22) 

This steep dependence on r, simply reflects the high suppression resulting from decoher- 
ence due to incoherent scatterings in the plasma whose frequency increases as 1/r. This 
suppression has the same origin as the one forbidding electroweak baryogenesis with CP 
violation originating from the mixing of light quarks in the MSM. For larger coherence 
time, Z-™ scales approximately as 

K oc ± (r-oo). (2.23) 

This fall-off with increasing r describes the semi-classical limit in which particles 
propagate a distance long compared to their Compton wavelength in which case, fast 
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oscillations of their wave-function wash away the interference required to generate a CP 
asymmetry. This fall-off takes place in the thick wall regime, r < w, the situation for 
which our derivative expansion applies; it persists until r ~ w, which defines the thin wall 
regime. The two behaviors ( |2.22| ) and (|2.23| ) are easily identified in Fig. 3. 

For the case of interest, the damping rate is essentially dominated by weak interaction 
processes: 7^ ~ ct w T . A crude estimate for the coherence time r~ h , ~ 7iT ' yields the 
range 20/T < r~ h < 30/T, which lies comfortably in the asymptotic domain described by 
(PI). 

In summary, an analytic form which fits well the source term for the Higgsino number 
in the domain of interest is Eq. (|2.20[) with 



22 1 e m / T /t^\ 5 / 2 

Xr « — — t j^r — 20/T < rr < 30/T . (2.24) 

h 5 Tr h (l + e m / T ) 2 \mJ I - h- I \ ) 

This form is only a fit which is valid in the range 0.5 < m/T < 1. 
Let us assess the domain of validity of the expression above. 

1. Results ( |2.21| )-( |2~2"3| ) are not to be trusted in the limits r > w where our derivative 
expansion does not apply. Instead, in this limit, one expects to observe the dependence 
on r to weaken as r — > w, and to vanish as r > w; that is, one expects to observe the 
factor T™ to saturate at a value r ~ w, to become a function of the mass only. 

2. The coherence time t$ given above is an estimate based on our knowledge on the 
damping rate which, so far, has only been studied in the low momentum limit p4f ; 
this value can only be a crude estimate. 

3. Finally, we turn to the important question of the validity of the mass expansion 
we used to derive formulas ( |2.20| )-( p.23| ). A careful study of this expansion shows 



that its expansion parameter is either mA in the limit wA < 1 or ra/w in the limit 
wA > 1. In both limits, the expansion parameter is less than one because of the rela- 
tion m < uj = \/k'j_+ m 2 . Let us give the physical interpretation of those statements!. 
In a typical scattering off a diffracting medium characterized by a step potential of 



For a more detailed discussion see Ref. 0. 
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height V, reflection and transmission amplitudes result from the constructive and de- 
structive interference of various diffracted waves generated everywhere inside the bulk 
of the medium. Transmission and reflection amplitudes will be comparable^ if the 
incoming wave penetrates coherently the diffracting medium over at least a distance 
of order 1/V, and has few oscillations over that distance. Suppression of the reflection 
amplitude arises if the coherence length A of the incoming wave is smaller than 1/V 
or if its energy uj is larger than V. In the case where A <C only a layer A of 

the medium effectively contributes to the coherent reconstruction of the reflection and 
transmission amplitudes, this is the phenomenon of decoherence, which suppresses the 
reflection amplitude with powers of A V. In the case where V C w, fast oscillations of 
the propagating wave inside the medium, tend to attenuate the reconstruction of the 
reflected amplitude with powers of V/u. In the present case, the "diffracting medium" 
is the wall and its height V is the mass m of the scattering particle, hence, suppression 
factors are controlled by A ro or m/w, whichever is smaller. An alternative method of 
computation of the currents J± consists of computing R and T by solving Majorana 
equations including the imaginary part of the thermal self-energy. This method au- 
tomatically accounts for both effects occuring here @. In particular, use of the mass 
expansion in this context suggests an expansion parameter ~ m A / yl + uj 2 A 2 , which 
corroborates the analysis above. 
- S quarks - 

We now turn to the calculation of the source for the axial stop number. The stop 
mass matrix is given in ( |2.7| ) and the top axial charge operator Qs is given in ( |2.9| ). As for 
the Higgsino number, we proceed in computing the current source J± in the wall frame 
using Eqs. ( |2.2|) , which we then input into formula ( |2.4|) to construct the source 75 (x,t). 

This time, the amplitudes are computed in solving a set of coupled Klein-Gordon 
equations. We obtain up to an overall phase, at leading order in M. 2 , 

10 In order to obtain significant CP violating asymmetries, both reflection and transmission 
amplitudes are to be significantly different from zero, otherwise, as either \R\ or \T\ goes to zero, 
the other one goes to one from unitarity, and both |i?| 2 — |i?| 2 and |T| 2 — |T| 2 vanish correspondingly. 
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We only displayed the contributions whose interference contribute to a CP asymmetry in 
J+. These specific paths are depicted in Fig. 4. As we did earlier, we assume the squark 
density matrices p Zo and p Zo +A to describe thermal distributions in the unbroken phase 

Dmg(n b AL (E, +v),n b AR (E, +v) 

Dmg^ L (E,-v),n b ^ R (E,-v) 

n h m (E, v) now, refers to the Bose-Einstein distribution, n&, boosted to the wall frame, 

n b m = (E X p[-f w (E - v w k±)] - l)" 1 . (2.28) 



Pz +A 



(2.27) 



The soft supersymmetry breaking masses are kept in ( 2.27|) as, for large values, they yield 
an exponential suppression of the baryon asymmetry produced. 
From ( |2.25| ), we obtain for the current sources 



J+ = (1, 0, 0,v) x | 



+ [ A dZl P dz 2 SklMz r Z2) xlmTr 
Jo Jo u 2 

f A , f A , sin 2u(zi — zo) T _ 
+ / dzi dz 2 ±± x Im Tr 



p Zo Qs Mj Ml 



p Zo +A M 2 2 Qs M\ 



and 



J- = J+(Pz <-> p Zo +A, v <-> -v). 
Performing an expansion in the wall velocity v w , we find, in first order in v z 



(2.29) 
(2.30) 



(J+ + J-T =o 

(J+ + J_)° = -y w v w x J s x ^ 



i=m L ,m R 



d 3 k g(u;A) 



(2tt) 3 4a; 5 
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(2^ 



r>Ei/T E . (2.31) 



(l_ e ^/T)2 T 



where g(£) is defined as 

0(0 = 1 - cos 2£ - £sin2£ 



(2.32) 



and Js is a new CP violating invariant given by 

Js = - d z a a 4 /T 5 

A 2 (2.33) 
= 4-|A/isuu> B -<j> A )M% v {T,z)d z p . 

After a few simple manipulations, we derive the following expression for the stop axial 
source jg(x,t), 



rp4 

+ 0{vl,(r/w) 2 ) 
where N c is the number of colors, = 3. The function I™ is given by 



If = ~ A ^T dyy 2 / dZ -) (2.35) 

it is simple to show that the factor If 1 vanishes rapidly with r 

I™ oc r 3 (r -> 0) , (2.36) 

as, in this limit, incoherent plasma scatterings become overwhelming. This behavior, 
already noted in the Higgsino case, is a universal property which can be traced to the 
quantum nature of CP violation, conflicting with the classical nature of the plasma physics. 
For larger coherence time, Xf 1 behaves approximately as0 

(r-oo). (2.37) 

Both behaviors (|2.36|) and ( |2.37|) are evident on Fig. 5. Unlike the case of the Higgsinos, 
the squark plasma physics is dominated by strong interactions, and so is the damping rate, 
hence we estimate t§, oc (2a s T) _1 , to be about 5/T. From this estimate, we infer that 



11 This behavior is cut off at a value r ~ L; at this value, Z|™ is expected to saturate to its 
"thin wall" value. 
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the regime of relevance is neither the decoherence regime nor the semi-classical regime 
but rather an intermediate regime corresponding to the peak shown on Fig. 5. This is a 
situation already encountered in Ref. ||, in the case of the top quark in the two Higgs 
model (cf. §4). 

In summary, and with the above value for Tg, an analytic expression which provides a 
reliable fit to the source 75 (x, t) for the axial stop number, in the range 0.5 < m/T , is0 

N C T 4 \ - T e m l T 
j s (x,t) ~ lwV w n x Js x > — -, 7777777 with T S ~ 5/T . (2.38) 

m=rht, jffiR 

This expression has been derived under the same assumptions as the ones made to derive 



the corresponding analytic form (|2.24|) for the source for the Higgsino number, 7^(x, t). 



These assumptions have been evaluated in the discussion following Eq. ( j2.24| ). 



3. Diffusion equations in the SSM 

Only those particle species which participate in particle number changing transitions 
which are fast compared with the relevant timescales, but which carry some charge which 
is approximately conserved in the symmetric phase, can have significant nonzero densities 
in the symmetric phase during the transition. If the system is near thermal equilibrium 
and the particles interact weakly, the particle densities rii satisfy 

m = kiHiT 2 /6 , (3.1) 

where \ii is a local chemical potential for particle species i, and ki is a statistical fac- 



tor defined by eq. (3J.). For light, weakly interacting particles ki ~ 2 (boson degrees of 
freedom) +l(fermion degrees of freedom), while for particles much heavier than T it is ex- 
ponentially small. If we consider a reaction which changes the particle number of particle 
species i by Aj, near thermal equilibrium the difference beween the rates for the reaction 
and its inverse will satisfy 

■p Si^iA^-p Ei n i^i fir /rp3 (n o\ 

J- net = 7p J- fluct = ^ t>-Lfluct/-t , {o.Z) 



12 We emphasize that this is only a fit. 
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where Tfl uct is the total rate for the reaction and its inverse per unit volume. For conve- 
nience we will henceforth define particle number changing rates to be (6/T 3 )Ffl uct . 

We can now write down a set of coupled differential equations which include the effects 
of diffusion, particle number changing reactions, and CP violating source terms, and solve 
them to find the various particle densities in the SSM. Anticipating a small departure 
from equilibrium, we incorporate particle number changing reactions and sources as two 
distinct terms. Diffusion is described by a standard diffusion term without a provision to 
account for the potentially fast relative motion of the sources in respect to the plasma. 
It is a good description in the regime of a wall velocity v w small compared to the speed 
of sound in the plasma c s = l/y/3. This condition, which is likely to be fulfilled in the 
minimal standard model |24|143[], may or may not be fulfilled in more general theories 
such as the ones considered here. To find out would require a complete calculation of the 
phase transition, which is beyond the scope of the present work. Further simplifications 
of these equations take place when we neglect all couplings except for gauge interactions, 



and the top quark Yukawa coupling. We include the effects of strong sphalerons [45 , 46 



but neglect the weak sphalerons until near the end of the calculation. The neglect of the 
weak sphalerons allows us to forget about leptons in our differential equations, and will 
turn out to be a good approximation when computing Higgs and quark densities. We also 
neglect the effects of hypercharge gauge forces and screening, which can be shown to affect 
the baryon number produced by a factor of at most order one [[IT|. The particle densities 
we need include q = (£2, + the right handed top quark t = tn, the Higgs particles 
h = {h~ + h° + h' + + h'°), and their superpartners q, t, h. The individual particle numbers 
of these species can change through the top quark Yukawa interaction, the top quark 
mass, the Higgs self interactions, and anomalous weak interactions, and the supergauge 
interactions. We will find that baryogenesis in the minimal model is only feasible if some 
of the superpartners of the gauge and Higgs bosons are light, so that we may take the 
supergauge interactions to be in thermal equilibrium (q/k q = q/kq,t/kt = t/k^,h/kh = 
h/kjj, and describe the system by densities Q = q + q, T = t + t and H = h + h. As 
shown in § 2, CP violating interactions with the phase boundary produce source terms 
7^ for the Higgsinos and j§ for the q — t densities, which tend to pull the system away 
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from equilibrium. When we include strong sphalerons (with a rate T ss ), we will generate a 
right handed bottom quark density, B = bn + bn, as well as first and second family quarks 
Q {1,2) Li > Ur , Cr , Sr , Dr. However since strong sphalerons are the only processes which 
generate significant numbers of first and second family quarks, and all quarks have nearly 
the same diffusion constant, we can constrain these densities algebraically in terms of B 
to satisfy 

Qil = Q2L = -2U R = -2D R = -2S R = -2C R = -2B = 2(Q + T) . (3.3) 

For simplicity we will also assume all squark partners of the light quarks are degenerate 
and take 

k QiL = k Q2L = 2k SR = 2k DR = 2k UR = 2k CR = 2k B . (3.4) 

We include scattering processes involving the top quark Yukawa coupling, with rate T y , 
and in the phase boundary and broken phase we have Higgs violating processes at a rate 
Th and axial top number violation at a rate T m . Following ref. ||, particle transport is 
treated by including a diffusion term. We take all the quarks and squarks to have the same 
diffusion constant D q and the Higgs and Higgsinos to have diffusion constant D^. 

The rates of change of the various densities are now described by the coupled equations: 

Q = D q V 2 Q - T y [Q/kq - H/k H - T/k T ] - T m [Q/k Q - T/k T ] 

- 6r ss [2Q/k Q - T/k T + 9(Q + T)/k B ] + 7s ~ 

f = D q V 2 T - T y [-Q/kQ + H/k H + T/hr] (3.5) 

- T m [-Q/kQ + T/k T ] + 3T SS [2Q/k Q - T/k T + 9(Q + T)/k B ] - Ts 
H = D h V 2 H - T y [-Q/kQ + T/k T + H/k H ] - T h H/k H + Th . 

Several simplifications of equations ( |3.5|) can be made. First we ignore the curvature of 
the bubble wall, and so T m , r\, and 7- ^ are only functions of z = |r + v w t\, where v w 
is the bubble wall velocity. We will assume that the density perturbations of interest are 
only functions of z, the coordinate normal to the wall surface. 

With these assumptions we arrive at the equations for Q(z), T(z), and H(z) in the 
rest frame of the bubble wall: 

20 



= -v w Q' + D q Q" - T y [Q/kq - H/k H - T/k T ] - T m [Q/kq - T/k T ] 

- 6T SS [2Q/k Q - T/k T + 9{Q + T)/k B ] + Is 

= -v w T' + D q T" + -T y [-Q/kq + H/k H + T/k T ] (3.6) 

- T m [-Q/kq + T/kr] + 3r ss [2Q/kq - T/k T + 9(Q + T)/k B ] - is 
= -v w H' + D h H" - T y [-Q/kq + T/k T + H/k H ] - T h H/k H + Th . 

We now assume that the rates T y and T ss are fast, and so Q/kq — H/kn — T/kr = 
0(l/T y ) 7 2Q/kq - T/k T + 9(Q + T)/k B = 0(1/T sa ). We will check later whether this 
assumption is self consistent. We then take the linear combination of eqs. ( |3.6| ) which is 
independent of T ss , T y , and substitute 



+0(l/T 3S ,l/T y ) . 



T = -H 



(3.7) 



k T (2k B + 9kg) 
k H ( k B + 9k Q + 9k T 

We then find that the Higgs density satisfies 

= -v w H' + DH" -TH + J + 0(1/T SS , 1/Ty) (3.8) 

where D is an effective diffusion constant, T is an effective decay constant and 7 is an 
effective source term, given by 

D q (9k Q k T - 2k Q k B - 2k B k T ) + D H k H (9k Q + 9k T + k B ) 



D 



9k Q k T - 2k Q k B - 2k B k T + k H (9k Q + 9k T + k B ) 



- / , w k H (9k Q + 9k T + k B ) . 

! " l '* ~ lh) \9k Q k T - 2k Q k B - 2k B k T + k H (9k Q + 9k T + k B ) > 



F=(r m + r 



h 



9k Q + 9k T + k B 



9k Q k T - 2k Q k B - 2k B k T + k H (9k Q + 9k T + k B ) 

In these equations, 7 is the sum of the rate of generation of axial quark number and 
Higgs number inside the wall as given in (|2.34j ) and ( |2.20| ) while, T is the total rate of 
relaxation for those charges. We estimate the latter to be 

(r „ + v h ) „ iM -^ A? S in^ + ■ (3-10) 
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Equation (|3.8|) is easily solved numerically for arbitrary shape of the source 7 and decay 
term F, however in order to qualitatively understand how the baryon number produced 
depends on the various parameters we will approximate the source as a step function of 
width w 

7=7, w > z > 

(3.11) 

7=0, z > w, z < , 
while for the decay terms we take 

f=f, z>0 

(3.12) 

r =0, z < . 

The effective diffusion constant is also spatially varying since the statistical factors ki 
depend on spatially varying particle masses and since the weak interaction cross sections 
depend on the Higgs vevs, however we will make the reasonable approximation that D 
is constant. An analytic solution to eq. ( |3.8| ), which satisfies the boundary conditions 
H(±oo) = is now readily found; for z < (the symmetric phase) this is 

with 

A = \ ^ = —T2 • ( 3 - 14 ) 

v w + J4DT + vt 



We will see that DT ^> and so a good approximation to eq. (|3.14|) is 



a*[1 ) U-e-^vi 



7s + T h 



(3.15) 



From the form of ( |3.13|) , we see that the CP violating densities are non zero for a time 
t = D/v^j, and so the assumptions about which rates are fast which were used to derive 
eq. ( |3.8|) are valid provided -Dr ss /v^, DTy/v^ ^> 1, DT^/v^, 1, and the scattering 
processes due to Yukawa couplings other than top are slow. 
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To estimate D we take the Higgs diffusion constant to be comparabie to the 
diffusion constant for feft handed leptons, which was estimated in the MSM in Ref. j|7 



to be 110/T and take D q from j|7| to be 6/T. (These numbers wiii decrease siightly due 
to the supersymmetric particie content of the plasma-we ignore this effect as being smaii 
compared with other uncertainties in our calculation.) For the k^s we assume that all 
the supersymmetric particles are heavy compared with T except for the neutralinos and 
charginos and so 

k Q « 6, k T ~ 3, k B « 3, kf{ m 12 . (3.16) 



We then find the effective diffusion constant defined in eq. ( |3.9|) is large, 

D « 100/T . (3.17) 

The large effective diffusion constant indicates that most of the transport of CP violating 
quantum numbers is done by weakly interacting particles, ie.the Higgs and Higgsinos, and 
since Yukawa interactions readily convert Higgs number into axial top number, transport 
of axial top number is surprisingly efficient. 

For the scattering rate due to the top quark Yukawa coupling we estimate 

« (27/2) A t 2 a s (C(3)Ar 2 ) 2 T (3.18) 

and so DTy/v^ > 2/t> 2 and the assumption that this rate is fast is self consistent. The 
next largest Yukawa coupling is the bottom quark's. Including scattering from this Yukawa 
coupling would give corrections to our results of order 

~ D ((27/2)A 2 a s (C(3)Ar 2 ) 2 T) /v 2 w . (3.19) 

We assume that the ratio tan f3 of Higgs expectation values is not unnaturally large, and 
so scattering due to the bottom and other Yukawa couplings may consistently be neglected 
for v w > 10" 2 . 

For the anomalous fermion number violating rates we take 

g 

T ws = 6ko4T , T ss = QK-a^T , (3.20) 
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where unknown parameters usually assumed to be of order one. Thus the weak 

sphaleron rate may safely be taken to be slow provided 

k/vI < 10 4 (3.21) 

and the strong sphaleron rate is fast if 

n'/vl > 5 . (3.22) 

In our computation of the baryon asymmetry we will approximate the strong sphaleron 
rate as fast and the weak sphaleron rate as slow. 

What we set out to compute was not the Higgs density in the symmetric phase but the 
total baryon number density left inside the bubble. We now turn the weak sphaleron rate 
on, assuming it has a negligible effect on particle densities (eq. ( |3.21| ) is valid), however 
it provides the only source for net baryon number. We thus take ps, the baryon number 
density, to be a function of z satisfying 

= D q p'g - v w p' B - Q(-z)n F T ws n L (z) , (3.23) 

where ul is the total number density of left handed weak doublet fermions, np = 3 is the 
number of families, and we have assumed that anomalous baryon number creation takes 
place only for z < (the symmetric phase). Eq. fl3.23|) has solution 



p B (- z ) = JJ^± [ nL ^)dz-^- I* dz'0(-z')n L (z')(l-e v ^*- z 'V D A , (3.24) 

Vw J-oo V w J V / 

which is a constant for z > and vanishes as z — > — oo. Thus, up to corrections of order 
r w8 D/v^, the baryon density inside the bubbles of broken phase is simply proportional to 
the integral of in the symmetric phase. 

We now return to eq. (|3.6| ), in order to compute n^. As pointed out by Giudice and 
Shaposhnikov ||46|| , if we use eq. (|3.16|) we will find the answer is zero in the limit T ss — > oo, 
so we need to compute the 0(1/T SS ) corrections to particle densities!!. We will assume 



13 In this limit, one should also include the contribution of sources for conserved charges: B — L, 
.... Local densities of conserved charges are also generated by the scattering of particles on the 
moving wall through charge separation. As discussed in section §2.1, these sources are subleading, 
however they do not suffer from the strong sphaleron suppression. 
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T y 3> T ss , (k' < 7) and take 



(3.25) 



Substituting these values into eq. fl3.6|) , we find 

o =D ,i<r + r») + «(Q' + to - 3 r„ g - L + ^) 

_ / -fc B (fc Q + 2/c T ) \ , „ _ , _ fkB±9kQ±9kr 
-\k H (9k Q + 9k T + k B )J [IJqli ~ VwU) 6lss { k B k Q 



(3.26) 



+ o(i/r ss ,i/r y ) 

D q H"-v w H'\ ( -k 2 B k Q (k Q + 2k T ) 



3k H (9k Q + 9k T + k B ) 



+ o(i/rJ„ i/r, 



We now solve algebraically for m = Q + Qil + Qil using eqs. (|3.3| ), ( |3.25|) , and ( |3.26|) . 
and find 

n L =5Q + 4T 

5k Q +4k T \ x ( 9k Q k T - %k B k T - 5k B k Q \ (3-27) 



/ 9/cq/;;t - - 5k B k Q \ 

Q V fc//(fcB + 9fc Q + 9fc T ) ; ' 



k Q / \ K H yK B -^yKQ 

If we use eq. ( |3.16|) we find 



n L = 75 Q = -1/56 I — — I 



(3.28) 



so the baryon density is proportional to r„, s /r ss , and is only sensitive to the ratio of 
k/k', provided eqs. ( p5.21| ) and ( p5.22| ) are satisfied. This reduces the uncertainty in the 
baryon asymmetry since estimates for both k and k' vary by several orders of magnitude 
but we expect the ratio to be approximately one. The result that ul and the baryon 
density are suppressed by a factor of l/r ss does not hold if one considers modifications to 
eq. ( |3.16| ) due to higher order corrections ]|6| or due to the contributions of nondegenerate 
squarks. The cancellation which makes the first term of eq. ( |3.27| ) dominate the second 
term no longer occurs when nondegenerate masses are considered. However examination 
of eqs. ( |3.26| ) and ( |3.27|) shows that corrections from this lack of cancellation are negligible 
unless either k'/v^ > 10 3 , or some squarks are not much heavier than T. Note that there 
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can be significant enhancement of the baryon density if, for example, only the top squark is 
light in the symmetric phase, as we will discuss at the end of this section. With all squarks 
heavy, our final answer for the baryon to entropy ratio in the broken phase, combining 
eqs. (CT) , ( |3TT5D , ( ggg ), ( ggfl ), and is 



Pb_ _ _ ( 3AT W8 \ ( 1 _D 1 

« - 3.5 x 10- 5 -f w v w ^) S 



(3.29) 



with 

S= (T fi^y^l . (3.30) 

We have taken the entropy s to be s = (2n 2 g*/45)T 3 = 55.1 T 3 (#* = 125 3/4) and 
we have made explicit the dependence on the velocity v w and the thickness of the wall w. 

The factor S is a dimensionless number, function of the supersymmetric parameters 
/U, sinc^B, wi,2, A, sin^, tan/3 and A/3, the total variation of f3 in the wall, as well as a 
function of the known gauge and top Yukawa couplings and W and Z masses. In short, 
S is a concise representation of the dependence of the baryon asymmetry produced on the 
yet unknown supersymmetric parameters of the SSM. 

To compute the baryon asymmetry, we need to compute the factor in parenthesis 
in Eq. (|3.30| ). This factor has its origin in the mechanism which transports the CP 
violating asymmetries in front of the wall. If this transport is efficient, the answer should 
become independent of the wall thickness w. Indeed, using our estimates, 100/T, for 



D given in ( |3.17[ ) and our estimates for V in ( |3.9[ ) and ( |3.10| ), we find 2yT/D ~ 5 x 
10 -3 a/7 sin 2 (3 + 1 T, a value fairly insensitive on the supersymmetric parameters. Hence, 
unless the wall thickness w is anomalously large,0 we find the factor in parenthesis to be 
equal to 5 x 10~ 3 a/7 sin 2 (3 + 1 and, at leading order, independent of the wall thickness. 

The iS-factor becomes, using general expressions for the source terms given in ( |2.20|) , 
( |2.21| ) and ( |2.34| ), ( |2.35| ), and our expressions for the relaxation rates Y m and Yh in ( |3.1(J| ), 



14 Typical estimates for wT range between 10 and 100. 

26 



S » 5.5 x 10" 3 ^ ^ sin 0s A^/ y7 sin 2 /3 + 1 x 

J £ ( 1 + ^Tn" ) + J f 2 + J f 1 ^ T7T ) sin <^ + ( 3 - 31 ) 
h m 2 10 h h m 2 10 / 

( [Jf - + If-] 10 A) si n(^ - d> A ) } 

The first term in parenthesis represents the contribution of the charginos and neu- 
traiinos while the second term represents the contribution of the top squarks. 

To try this formula, we use light neutralinos and charginos: (x = rhi = 2rh\ = 50 
GeV, relatively heavy squarks; ttil,r = 150 GeV, and take A = 50 GeV. We need to 
know the ratio of these masses over the transition temperature T. The latter has its value 
completely determined by the parameters of the theory; in our analysis, however, it is a 
free parameter. As an indicative value, we choose 60 GeV 0. For these values, we find 
It = Tf 2 « 5.5 x 10" 2 , J™ 1 « 3.5 x 10" 2 and lf L ' R « 2.4 x 10" 4 . the 5-factor becomes 

h h h s 

S « (-) ; Af3 (0.032 sin6 B + 5 x 10" 3 sin(0 s - <b A ) ) 5.5 x 10" 3 

W v/l + 7sin 2 /3 

« (1.8 x 10" 4 sin0 B + 6.6 x 10" 6 sin(0 B - A ) ) A/3 . 

(3.32) 

It is clear from the above equation that the squark contribution is only significant in the 
limit <pA ^ 4>b or in the limit the charginos and neutralinos are heavy. 

Gathering all the above information, we find the following results. The largest contri- 
bution arises from light charginos and/or neutralinos, in which case, the asymmetry can 
be as large as 



- J w v w (4) sin^s A /3 6 - 5 x 10 ~ 9 - ( 3 - 33 ) 



PB „ ( K \ -j, A a a k « m-9 

S V/t' 

The measured baryon asymmetry is (4 — 7) x 10 -11 . So, electroweak baryogenesis is 
significant provided that 



(^) \lwv w sin 0s A/3| > 7.5 x 10" 3 (3.34) 



15 Generically, one expects the temperature to be below the one in the SM ( ~ 80 — 100 
GeV), as, in the MSM, the superpartners contribute to the effective potential in a manner which 
decreases the critical temperature. 
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and is negative in sign. Let us discuss the magnitude of each term separately, k and 
k' are two not well-known parameters characterizing the strength of the electroweak and 
strong anomalous processes, respectively, however, their ratio is expected to be of order 
one. In the minimal standard model, the wall velocity, "f w v w , is no smaller than 0.020 
and is more likely of order 0.10 or larger [p4| , |43[| ; although no calculation has been done 
for SSM, it is reasonable to assume similar values. Finally, A/3 is the overall variation of 
the ratio of the two Higgs expectation values i>2 and v±. As we argue in §5, its presence 
is an artefact of working at fourth order in the mass, it can be removed at the cost of 
introducing additional mass suppressions. From these considerations, we infer an optimal 
bound on the CP violating phase <pB 

|sin0 B | > 0.025 . (3.35) 

Only with this bound satisfied, is electroweak baryogenesis achievable in SSM with light 
charginos and/or neutralinos and heavy and degenerate squarks. 

In the case of neutralinos and charginos which are heavier than T or </>^ ^> (ps ~ 0, 
only top squarks contribute to the asymmetry. If all squarks are degenerate, they must all 
be heavier than ~ 150 GeV in which case, the requirement becomes 

[A/3 sin(0 A -0 B )| > 0.65. (3.36) 

which leads to experimentally ruled out electric dipole moments. Also, unlike the 
charginos/neutralinos case, A/3 must be non-zero and its presence is not an artefact of our 
approximations; A/3 can be significantly smaller than one. Clearly, electroweak baryogen- 
esis in SSM with all neutralinos, charginos and squarks heavy is likely to be incompatible 



16 This lower bound corresponds to the situation of maximal damping of the motion of the wall 
in the plasma, that is, it corresponds to the thin wall situation where mean free paths are larger 
than the thickness of the wall, w. This lower bound is a decreasing function of the Higgs and 
top quark masses; the specific value 0.02 has been computed following Ref. [24 for the values 
run ~ 65 GeV and rat ~ 175 GeV. 

17 This larger value accounts for thermal scattering within the wall as, in the MSM, the wall 
thickness w is typically larger than the mean free paths of the W's and Z's, t w , and of the top 
quark, r s . Large uncertainties arise from our imprecise knowledge of the ratios r WtS /w. 
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with constraints from electric dipole moments. We will discuss bounds ( |3.35| ) and ( |3.36| ) 
along with their uncertainties in the last section. 

These conclusions are altered considerably if, say, the left handed bottom squark and 
left and right handed top squark masses squared are rather light, but the other squark 
masses are heavy. (This mass pattern could be a result of renormalization due to the large 
top Yukawa coupling near the Planck scale.) Then the factor multiplying H in eq. ( |3.27| ) 
does not vanish. If we take 



we have 



k Q « 18, k T « 9, ks ~ 3, ku ~ 12 , (3.37) 
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n L = , (3.38) 
D « 72/T , (3.39) 



and 



p B _ ( 81 ADT t 

s 



(3.40) 



V 82v 2 w s 

i.e., pb is enhanced by a factor of ~ ld>DT as /v^ over the case with no light squarks and 
fast strong sphalerons, and is sensitive to the weak sphaleron rate rather than the ratio of 
weak and strong sphaleron rates. After a few substitutions, we obtain 

P — « -1.5 x 10" 4 S/v w , (3.41) 
s 

where we have made use of the 5-factor defined in ( |3.30| ) and ( |3.31| ). To obtain a numerical 
estimate, let us assume the values A ~ p ~ rri2 — 2toi ~ 50 GeV, and ~ T ~ 60 

GeV. The factors in flOg ), are now 1? = if 2 w 5.5 x 10" 2 , if 1 w 3.5 x 10" 2 and 
X™ L - H w 0.018. the 5-factor becomes 

5 ^(-) 2 A/3 (0.032 sin0 B + 0.093 sin(0 s - <f> A ) ) 5.5 x 10" 3 

6 v / l + 7sin 2 /? . (3.42) 

w ( 1.8 x 10" 4 sin 0s + 5 x 10" 4 sm(<p B ~ <1>a) ) A/3 

This time the contribution of the top squarks - the second term in parenthesis, is potentially 
as significant as the one of the light charginos or neutralinos. Combining this result with 
(CT) , we find 
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— « (A/3 sin0 B 2.5 x 10" 8 + A/3 smU B - <f> A ) 7.5 x 10" 8 ) . (3.43) 

s v w 

This is a significant contribution to the baryon-to-entropy ratio provided that 

| A/3 sin 0s + 3 A/3 sin(0B - A )| > — 1.9 x 10" 3 . (3.44) 

and is negative in sign. 

If we use the range 0.1 — 0.3 for the wall velocity v w and the range 0.1 — 1 for k, 
and take masses for the superpartners which are optimal for baryogenesis, we obtain the 
following constraint on the magnitude of the phases 4>a and b 

light charginos/neutralinos/top squarks |A/3 sin 0b | > 2 x 10 -4 — 6 x 10 -3 
light top squarks/charged higgs | A/3 sin (0 b — 0a)| > 7x 10~ 5 — 2 x 10~ 3 

(3.45) 

We emphasize that in the term contributed by the charginos and neutralinos, A/3 will not 
be present at higher order in the mass expansion. By taking the top squarks to be light 
we obtain a possible two-order of magnitude enhancement in p B /s over the situation with 
all squarks degenerate and heavier that T (cf. ( |3.34|) - (|3.36|) ). 

We also obtain qualitatively similar results to eq. ( |3.40| ) if strong sphalerons are slow, 
i.e. eq. ( |3.22|) does not hold. Our formulae are also radically modified if weak sphalerons 
are sufficiently fast and/or if the wall velocities are so slow that eq. Q3.21|) is violated. 
Then most of our simplifications of the rate equations, such as the neglect of leptons, are 
invalid. We then expect the final answer for ps to be insensitive to the sphaleron rates, 
being determined by near-equilibrium physics. 

4. Baryon density in the Two Higgs model 

We can now easily solve for the baryon density in the two Higgs model since the 
particle transport equations are very similar to those in the SSM. Eq. (|3.6| ) is unchanged, 
if we take the squark and Higgsino contributions to be zero, 75 to be the source for axial 
top number due to the top quark, 7 9 and we substitute 7b, the source for Higgs number 
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due to the Higgs particles, for 7^. Finally, in the case of two (one) light Higgs, the statistical 
factors become 

k Q « 6, k T ~k B ^ 3, k H « 8(4) , (4.1) 
and the effective diffusion constant, from ( |3.9|) 



^ 96 /88\ 

In Ref. § we computed 7 9 to be 

lq {x,t) ~ -— | Tut ; w T |m t | 2 z + («2 , (r/«;) 2 ) , (4.3) 

where m t (z), = \m t (z) \ e ld ( z \ is the space-dependent mass of the top quark expressed in 
the wall frame. To find 7#, we need to track the evolution of the Higgs number carried by 
the Higgses Hi and Hi as they evolve in the background of the wall. The space-dependent 
mass matrix is, in the basis Hi, H±, H% and H 2 , 



/ ... m^ 1 e~ i011 ... m\ x e~ id ^ \ 
miie* 1 . . . mj 2 e iei2 

m\ 2 e~ idl2 . . . ml 2 e~ i02 



(4.4) 



Vm^e^ 21 ... m\ 2 e id22 ... / 

We only displayed entries which violate Higgs number as they are the ones which control 
the charge generation as Higgs particles flow across the wall. The Higgs number charge 
operator takes the form 

Q H = Diag( 1,-1, 1,-1). (4.5) 

The analysis follows the steps of the one of the stop axial charge generation; in particular, 
Eq. Q2.29f) and fl2.30|) are directly transposable. We obtain 



(4.7) 



lH (x,t) = lw v w JL x {jp T£ + jp I]!} 2 ) + 0(v 2 w ,(t/w) 2 ) , (4.6) 
with 

Jg 1 = 2d z 8n mfi/T + d z 9 2 i rr^JT + d s B 12 m\ 2 /T 
Jp = 2d z 9 22 m\ 2 /T + d z 02i mii/T + <9 Z 12 m\ 2 /T . 
The function is identical to the one computed for the squark, I- 1 in Eq. ( |2.35|) . The 
damping rate is set by weak interactions, our estimate is tjj ~ 25/T. We choose for the 
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on-shell masses mi and m 2 of the propagating Higgses the zero-momentum contributions 
that the Higgs particles receive from plasma interactions in both phases: m\ ~ m 2 ~ T/3 

With these values, X^ 1 ph 0.25 and, within our approximations, 



lH {x,t) ~ lw v w ^ X 

{d z 9 n mjjT + d z 9 22 mt 2 /T + d z 6 21 m 4 21 /T + d z 9 12 m\ 2 /T) . 



Combining eqs. (ETI3D , ( gTj) , ( ggj ), (|3~26l) , and fl3~2"7D , the baryon to entropy ratio 



is 

p_b Z' 3^.r ws "\ ( d 



s v 5 6sr ss y v D j ^ 4 g^l 

« - 2.4 x 10- 5 7 ^ w (^) ft 
with most of the parameter dependence contained in the Ji factor 

(r m + I^^i^T 2 1 wT 



From (|3.9|) , we compute f = 0.11 (T m + Th). We estimate r m « A 2 T/21 and we 
parametrize rs A 2 T/140 with A, an undefined parameter function of the Higgs quartic 
couplings. As in the SSM, charges diffuse a long distance in front of the wall D ^> T and 
the term in parenthesis in Eq. ( |4.10| ), is largely independent of the wall thickness w. We 
find 

H = (7g + 7fe) w 7 x ID" 3 . (4.11) 

y/(T m + T h )/T^ w v w T^ 

Without going into a difficult study of the vast parameter space of the two Higgs models, 
we can obtain a fair estimate of the above quantity by neglecting the Higgs contributions 
to the source and to the rate, for the following reasons. The source 'jn written in (|4.8|) , 
is a linear combination of terms of the form m\^lTd z Qij. These terms all violate Higgs 
number, hence, they are proportional to the quartic self-couplings ~ X 2 (A9ij/w)T and 
are smaller than the contribution from the top quark ~ \ 2 (\A9ij/w)T 4 , unless the Higgs 
sector is strongly coupled. Similarly, we expect T m > Th- Under these assumptions, 

H » -10" 2 A6 (4.12) 
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and 



Pb 



s 



w 



( 



— ) A9 2.5 x 10 



k'J 



-7 



(4.13) 



Choosing the illustrative value 'fwVw ~ 0.3, this baryon per entropy ratio is significant 
provided that 



5. Outlook 

5.1. Accuracy of present computations of the baryon asymmetry 

We now look back on the many approximations and uncertainties present in our 



- Approximations - 

1. For the purpose of solving Majorana, Dirac and Klein-Gordon equations, we performed 
an expansion in powers of A4(x, t). The benefit was to work analytically and to 
express the answer as a sum of CP violating invariants. The convergence of this 
expansion has been discussed and established in the discussion following Eq. ( |2.24j) . 
We have further approximated the density matrices describing particle distributions 
with Bose-Einstein and Fermi-Dirac distributions for on-mass shell particles in the 
unbroken phase. In ignoring the non-equilibrium component of the distribution, we 
are ignoring corrections of order v 2 ^ (cf. discussion following Eq. (|2.13| )). In assuming 
on-mass shell particles in the unbroken phase, we are ignoring corrections of order 
(m/T) 2 and (gv/T) 2 . Some particles, such as the squarks, are expected to have 
SU(2) x £7(1) symmetric contributions to their masses which may be larger than 
the critical temperature (T c ~ 50-100 GeV). However, as discussed in §5.2, heavy 
particles do not contribute significantly to baryogenesis. So, at best, we expect that 
accounting for the full mass dependence yields numerical corrections of order one. 
One exception is that for the neutralino and chargino contribution, when we work 
to lowest nontrivial order in the masses we obtain a result proportional to A/3-the 
change during the transition in the angle specifying the ratio of the Higgs vevs. There 




(4.14) 



analysis. 
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is no reason to expect this suppression factor to persist at higher orders in a mass 
expansion. 

2. We defined our sources 7q in a layer of a size r, the coherence time. To postpone 
recourse to numerical methods, we assumed r to be smaller than the wall thickness 
w and ignored corrections of order (t/w) 2 . This is a very good approximation for 
strongly interacting particles but not necessarily for weakly interacting particles for 
which r is in the range (20 — 30) jT while the wall thickness w can span the interval 
(10 — 100)/T. Only a precise calculation of these two quantities can decide the quality 
of this approximation. The largest (t/w) 2 corrections are contained in the factors 

-. As explained in § 2, we expect these factors, which at most increase linearly 
with r, to "saturate" for r ^> u>, at about their values at r ~ w. For this reason, we 
do not expect higher order terms to bring large corrections to our analysis. 

3. We have made a number of simplifications of the equations describing particle trans- 
port and number changing processes. First, we assumed that deviations from thermal 
equilibrium were sufficiently small to allow us to describe particle distributions in 
terms of local chemical potentials and to make a diffusion approximation to transport 
processes. We expect this assumption to be quite good in the weakly interacting mod- 
els considered. We simplified our treatment of diffusion in neglecting the finiteness 
of the speed of sound, that is, we worked at leading order in an expansion in v w /c s 
where, c s = Our choice of the magnitude of the wall velocity is such that it is 
a fair approximation. Should the wall velocity approach or be larger than the speed 
of sound, diffusion is not a good approximation to transport and our computations 
are invalid. An improved calculation which covers large wall velocities has yet to be 
developed. We made a severe approximation by simplifying the wall shape (eqs. ( |3.11| ) 
and (|3.12|) ), which we expect to give an 0(1) estimate of the true solution. We also 
made assumptions about the approximate rates of strong and weak sphalerons and the 
wall velocity, (eqs. ( |3.21|) and ( |3.22|) ), and simplified our equations by assuming that 
the interactions proportional to the top Yukawa coupling were in thermal equilibrium. 
The size of corrections from these assumptions depends on how well the inequalities 
( |3.21| ) and ( |3.22| ) are satisfied. We also assumed similar diffusion constants for all 
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quarks, an error of order o^/ct^ ~ 10%. In fact this approximation for the quark 
diffusion constants is of very small numerical significance since diffusion is actually 
dominated by the weakly interacting Higgs, which provides a local source for axial 
top number far from the bubble wall. We also gave approximate estimates for the 
statistical factors k{ defined by eq. (lOp-here we expect corrections of order a few 
percent for light particles. The corrections to the quark statistical factors are impor- 
tant if k'/vI > C(10 2 ) since from eqs. ( |3.24|) , Q3.27|) and ( |3.28| ) we see that they give 



the only contribution to the baryon asymmetry which is not suppressed by the strong 
sphaleron rate. 

4. In most cases, our CP violating particle sources are dominated by the large CP vio- 
lation in the transmission of low momentum particles over a distance A whose wave- 
length is comparable to A. For these particles, kinetic theory starts to break down, 
giving corrections of 0(1) to our treatment. 

5. We have neglected the effects of long range gauge fields, which in general have an 0(1) 
effect on the baryon density |39| - |41[ . 



6. We have not included the contributions of the transport of conserved charges (such 
as B-L) to the baryon asymmetry. Such effects are higher order in wall velocity and 
masses, but may not suffer from the strong sphaleron suppression. We expect inclusion 
of such effects to change our results by at most 0(1). 

- Uncertainties - 

Uncertainties in our estimate of the baryon asymmetry reflect not only the approxima- 
tions above but also, and dominantly, our poor knowledge of certain parameters. Those are 
the coherence times r- ^ and t#, 9 , the diffusion constants D, the reaction rates V y ^,m and 
the parameters k and k' measuring the strength of the anomalous processes. Fortunately 
in the most interesting situations, the latter occur in ratio which significantly decreases the 
uncertainty in the baryon asymmetry. Also, for large k the baryon asymmetry becomes in- 
sensitive to n. Much work is needed to refine the determination of these parameters. Only 
r, k and k' require understanding new physics; the determination of the other parameters 
faces only technical challenges. The parameters v w and w describing the phase transition 
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are also left free, both because they are parameter dependent and because of the lack of 
accurate computations for the models under consideration. 

Finally we come to the main uncertainty, which is our lack of knowledge of the correct 
model of weak symmetry breaking and of the many new parameters introduced by any 
extension of the MSM. It is our hope that computation of the baryon asymmetry can 
provide a useful constraint on the weak symmetry breaking sector and on CP-violation. 

Because of the above uncertainties and also because of the approximations that we 
described earlier, we believe that the computability of the baryon asymmetry produced is 
reliable to an order of magnitude. It is with this caveat that we now present our conclusions. 

5.2. Can the baryon asymmetry be produced in the SSM? 

Previous work on baryogenesis in supersymmetric models neglected the enhancing 
effects of transport, and concluded that sufficient CP- violation for baryogenesis in super- 
symmetric models could be marginally consistent with electric dipole moment constraints 
if one made optimistic assumptions about baryon number violating rates in the phase 
boundary HH, and if chargino and neutralino masses were not too heavy. Our work 
shows that with reasonable assumptions about the rates of anomalous processes, sufficient 
baryon asymmetry can be produced with small CP violating phases of order lO - ^ -4 **, pro- 
vided that the top squarks and either the neutralinos or the charginos are light compared 
with the transition temperature. If only the top squarks are light, it is also required that 
the ratio of Higgs vevs is not fixed during the transition, while when the inos are light, 
the lowest order contribution in m/T is suppressed unless the ratio of Higgs vevs changes 
during the transition. The latter requirement implies that the effective theory during the 
transition has more than one light Higgs, which in turn means that at zero temperature the 
pseudoscalar and charged Higgs masses are not extremely heavy compared with the lightest 
Higgs mass. A light charged Higgs makes a potentially ruled out contribution to b — > 57 
||48|| unless partially cancelled by a contribution from a loop containing light charginos and 
stops. We conclude that as far as sufficient CP asymmetry is concerned, the SSM with 
some light superpartners ( < 100 GeV) is a good candidate for baryogenesis. With light 
superpartners and with CP-violating phases of order lO - ^ 2-4 -*, neutron and atomic electric 
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dipole moments will be below the current experimental bounds [33, 35 1 . Furthermore, a 
large fraction of the relevant range of masses for the superpartners coincides with the range 
to be probed by LEP II. 

In a calculation assuming only one light Higgs, the MSSM, with minimal superpartner 
content, has been shown to produce a phase transition sufficiently strongly first order to 
preserve the baryon asymmetry only when the lightest Higgs boson mass is less than 70 



GeV and when at least some of top squarks are lighter than 110 GeV ||36|| . We do not 
expect these bounds to be weakened significantly when the full parameter space for the 
Higgs masses is considered. Thus the baryon number washout constraint on the MSSM 
seems more powerful than the constraint of sufficient CP-violation. It is however subject 
to the uncertainties in the perturbative calculations of phase transition parameters. 

5.3. Conclusions about the baryon asymmetry in the two Higgs model, and comparison 
with other calculations 

To summarize § 4, sufficient baryon asymmetry may easily be produced in a general 
model with two Higgs doublets and soft CP-violation in the scalar potential, with CP- 
violating phase as small as ~ 7 x 10 -4 (eq. ( |4.14| )). This result allows for a much smaller 
phase than most earlier calculations in the two Higgs model. Here we explain how our 
calculation differs the earlier ones. 

The baryogenesis mechanism of axially asymmetric top quark reflection from the bub- 
ble walls |16[] also allowed a small phase of order 10 -5 in the two Higgs model, but only 



for the fine-tuned case where the bubble walls were thin, of order the inverse top mass. 



Refs. [|16| , |31fl concluded that when the walls are thick, a completely negligible CP-violating 
asymmetry is produced in the symmetric phase from top quark reflection. However in 
those papers several significant effects are neglected, such as thermal scattering within the 
phase boundary which is especially important for thick walls. Thermal scattering pro- 
cesses tend to interfere with baryogenesis by destroying the quantum coherence necessary 
for CP- violation []7|U, but also can in some cases enhance the baryon asymmetry produced. 
The enhancement comes about because CP-violating charge expectation values within the 
bubble wall can be converted to CP- violating thermal particle distributions inside the wall 
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by incoherent thermal scattering processes, and these CP-violating thermal particle dis- 
tributions can then diffuse into the symmetric phase, where they bias the relatively rapid 
anomalous weak processes towards producing net baryon number. We therefore find that 



the huge suppression of the top quark contribution to baryogenesis, found in refs. [ |16| , |3~1 



when the bubble walls are thick, is absent when thermal scattering and transport processes 
are considered. Instead, we find that the baryon asymmetry is not very sensitive to the 
width of the boundary. 

Let us now compare our method of computation to two alternative methods which 
have appeared in the literature. For the case of thick boundaries an alternative method of 
calculation of the particle distributions in the wall, which should be about as accurate as the 



thick wall approximations we made, would be to use the method of linear response [18,49 



i.e. to compute the charge current density produced from an initial CP symmetric thermal 
particle distribution when space-time dependent CP violating terms in the Hamiltonian 
are turned on for a time equal to the thermalization time r, and then dividing by r to get 
the rate for production of a CP-violating charge in the phase boundary. Such a calculation 
can be done diagrammatically, e.g. by computing the diagrams considered in ref [|5(| (which 
however does not contain a linear response calculation, as in that work the effects of a finite 
r are neglected). If one considers times longer than r in the linear response, including 
the damping terms in the quark propagators which are generated by gluon exchange is 
essential. Another method of calculation has been developed in Ref. [SI]. It consists of 
writing a Boltzmann equation for a one-particle distribution function which incorporates a 
CP-violating force term arising from the CP-violating space-dependent background. This 
equation is then solved for the resulting departure from thermal particle distributions, 
which is asymmetric between particles and their CP conjugates, and which are described 
in terms of chemical potentials. The latter are, in turn, inserted in a rate equation to 
compute a baryon asymmetry. It is not clear to us how to generalize this method to cover 
the case where several species mix, as occurs with the neutralinos and higgsinos in the 
SSM. In any case this method is semiclassical in nature and only describes the regime 
of r much larger than particle wavelengths (Tr ^> 1). In this regime our calculations 
also produce a semi-classical fall-off of the CP violating sources (cf. Fig. 3 and Fig. 5), 
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in qualitative agreement with the analysis of ref ||51|| ll3. This calculation is appropriate 
for weakly interacting particles such as the r-lepton in the two Higgs model, but not 
for strongly interacting particles whose mean free paths are not much longer than their 
wavelengths and for which decoherence effects are already perceptible. 

It is now evident that transport processes, omitted in the original thick wall calcu- 
lations P,|T7]rn|, significantly enhance the baryon asymmetry produced during the weak 
transition. In fact it has been suggested that the r-lepton plays a leading role in baryo- 
genesis due to its large diffusion constant P^| , |31| . However axial top quark number is also 
efficiently transported, because the large top Yukawa coupling allows axial top number to 
convert to Higgs number, which is transported by weakly interacting Higgs particles. An- 
other argument in favor of the tau lepton contribution to baryogenesis dominating that of 
the top quark is that the the axial top number tends to be washed out by strong sphaleron 
processes. In fact we find that this suppression factor is only about ~ 1/50 for k' of order 
one and wall velocities v w ~ 0.3. Furthermore, even for arbitrarily fast strong sphaleron 
rate, the strong sphaleron suppression will never be more than about 10 -3 , due to the non- 
degenerate thermal masses of the quarks |[46|| . Despite the suppression factors for the top 
quark contribution, we believe the tau is likely to be less important than the top for baryo- 
genesis in two Higgs models, because the source for axial tau number is suppressed relative 
to the axial top source by a factor of A^/A^, which is about 10 -4 unless tan/3 is large. 
In the SSM, it is only possible to avoid having sphalerons wash out the baryon number if 
tan P is relatively small f3(| , and so there is no significant effect from the tau or scalar tau. 
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These authors coined the word "classical baryogenesis" to describe their analysis. However 
their "classical force" is not completely classical. That is, the spin dependent CP-violating term 
in their one particle Hamiltonian is oc hd z 8, where 9 is the argument of the top quark mass. 
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Figure Captions 



Fig. 1. (a) Amplitudes contributing to J + . (b) Amplitudes contributing to ,/_. 

Fig. 2. (a) Contributions, to order (.M/T) 4 , to the transmission amplitude T of the neu- 
tralinos and the charginos. (b) Corresponding contributions, to order (.M/T) 3 , 
to the reflection amplitude R. 

Fig. 3. (a) The factor T™ plotted versus tT. Z~™ contains kinematic information on the 
propagation of the neutralinos and charginos in the plasma, (b) Its dependence 
on the mass eigenvalue m. The dots are the result from numerical integration 
and the solid lines are the fit (]2.24|). 



Fig. 4. (a) Selected contributions, to order (.M/T) 4 , to the transmission amplitude T 
of the squarks. (b) Leading contributions, to order (.M/T) 2 , to their reflection 
amplitude R. 

Fig. 5. (a) The factor If 1 plotted versus tT. This factor contains kinematic information 
on the propagation of the squarks in the SSM and on the propagation of the 
Higgs particles in two Higgs models, (b) Its dependence on the mass eigenvalue 
m. The dots are the result from numerical integration and the solid line is the fit 
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